
import numpy as np
#from matplotlib import pyplot as plt
#from scipy import ndimage
import sys 
nx=540
ny=540
nz=540
model=np.ones((nx,ny,nz),dtype='float32')
model=model*2000;
#model[:,:,0:100]=1000;
#model[:,:,100:230]=1500;
model.tofile("./model/model.dat_1")
#plt.imshow(model.T)
#plt.show()
